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PROPULSION SYSTEM MATHEMATICAL MODEL 
FOR A LIFT/CRUISE FAN V/STOL AIRCRAFT 


by 

Gary L. Cole and James F. Sellers (Lewis Research Center) 
and Bruce E. Tinling (Ames Research Center) 

SUMMARY 


A propulsion system mathematical model is documented that allows 
calculation of internal engine parameters during transient operation. 
A non-real-time digital computer simulation of the model is also 
presented. It can be used to investigate thrust response and 
modulation requirements as well as the impact of duty cycle on engine 
life and design criteria. Comparison of simulation results with 
steady-state cycle deck calculations showed good agreement. The model 
was developed for a specific 3-fan subsonic V/STOL aircraft 
application but it could be adapted for use with any similar 
lift/cruise V/STOL configuration. 
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INTRODUCTION 


A V/STOL aircraft propulsion system must supply the necessary lift 
forces and control moments during hover and vertical operations. To 
provide the desired aircraft handling qualities the propulsion system 
necessarily becomes an element in one or more high gain control loops. 
The result is a . potential for undesired interactions between the 
flight control and the propulsion system. In the case of multiple 
engine systems there must be a means of power management and, 
possibly, accommodation of an engine failure. 

In order to satisfy the unique requirements of V/STOL the propulsion 
system will have to meet certain thrust response, modulation and 
precision setting specifications. The harsh duty cycle of a V/STOL 
and the associated manipulation of the propulsion system will have an 
impact on engine life. 

Mathematical models are certain to play an important role in the 
investigation of V/STOL propulsion requirements. A simple 
linear-transfer-function propulsion model is not adequate. Although a 
detailed aero thermodynamic model is not necessary it is desirable to 
be able to examine some internal engine parameters (e.g. turbine inlet 
temperature, compressor exit pressure). 

The propulsion system mathematical model presented in this report was 
developed as part of a joint program with NASA Ames to model the 
research and technology (RTA) V/STOL aircraft shown in figure 1. The 
objective of the program was to provide data whictT could be used to 
investigate propulsion-system/flight-control interactions and pro- 
pulsion requirements for subsonic V/STOL aircraft. This was accom- 
plished by simulating the approach trajectory up to or just before 
hover while the aircraft was under automatic control. Approaches 
were made in the presence of disturbances such as turbulence, ini- 
tial vertical and lateral offsets and engine failures. A non-real- 
time digital computer simulation of the models was used. The air- 
craft and flight control system models, some flight path results, 
and availability of data are discussed in reference 1. Additional 
information regarding the RTA characteristics and modeling are 
given in reference 2. 

The objective of this report is to document the details of the 
propulsion system mathematical model as it was used in the RTA 
simulation. Although the model was developed for the RTA application 
it could be adapted for use with any similar lift/cruise "V/STOL 
configuration. The model is nonlinear and allows monitoring of 
internal engine pressures and temperatures. Dynamic representations 


1 



of rotor inertias, heat soak, fuel control and pitch actuators are 
included in the model. Steady-state accuracy of the computer 
simulation is discussed and some transient resultB are prosonted. 
Listings of the digital computer program are also given with a brief 
discussion of each subroutine's function. 


PROPULSION SYSTEM DESCRIPTION 

r 


The research and technology aircraft propulsion system consists of two 
lift/cruise turbofan engines, one turboshaft engine and one remote 
lift fan as shown in figures 1 and 2. The core engines are modeled 
after modified Detroit Diesel Allison XT701-AD-700 engines and the 
fans are based on Hamilton Standard 157.5-cm diameter variable pitch, 
low-pressure ratio fans. Additional information regarding the design 
of these units is given in references 3, 4, and 5. The lift/cruise fans 
are driven by the corresponding turboshaft engine low pressure (power) 
turbine through a reduction gear assembly. All three fans are 
connected by shafting through a combiner gearbox. This allows power 
transfer and prevents a loss of fan operation and hence thrust in the 
event of an engine failure. When a failure occurs, the failed engine 
power turbine is disconnected by means of an overrunning clutch to 
minimize the power loss. The remote fan is disengaged by declutching 
during conventional flight. Thrust amplitude is modulated primarily 
by varying fan-blade pitch angle. Total uninstalled thrust with all 
engines and fans operating at intermediate power is about 165 KN. 
Thrust is vectored by means of hooded nozzles on the aft engines and a 
louver system on the remote fan exit (fig. 1). The vectoring was 
included as part of the RTA airframe model. 


PROPULSION SYSTEM MODEL DESCRIPTION 


Engine Component Models 


This section describes some of the propulsion mathematical model 
details and assumptions and, where necessary, methods of 
implementation. 

Listings of the computer simulation are given in appendix B along with 
a brief description of each subroutine. The model equations can be 
readily determined from the FORTRAN program and are not summarized 
elsewhere. Definitions of the FORTRAN variables are given in 
appendix A. 
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A schematic of one turbofan engine is given in figure 3 and shows 
engine station numbers. The same station numbering is retained for 
the separate turboshaft and remote fan units unless noted otherwise. 
A computational flow diagram of the propulsion system model is given 
in figure 4. Each propulsion unit has its own representation. 
However, only two representations are shown in figure 4 because the 
turbofans, units 1 and 2 (fig* 2), are identical. 

Required aircraft inputs to the model are altitude and Mach number. 
Altitude is used as the independent variable to determine ambient 
conditions from standard atmosphere tables. An additional input, DTT, 
is added to ambient temperature to simulate nonstandard-day 
temperature conditions. Mach number is used to compute free-stream 
total temperature and pressure. In this model the inlet is simply 
treated as having constant pressure recovery. Ordinarily recovery 
will depend on air speed and angle of attack (see ref. 6 for example). 
This is a more important consideration in a tilt nacelle application, 
especially in the flight regime where air speed is high (120 Kts) and 
angle of attack is high (60°). 

Fan temperature and pressure ratios and corrected airflow are found as 
functions of both fan blade pitch angle (ranging from -^20° to 7.3°) 
and fan corrected speed (ranging from 70% to 110% of design) . The 
lift/cruise fan pressure and temperature ratios are different at the 
hub (station 25) and tip (station 13). A fan stall-margin calculation 
is made based on a knowledge of fan corrected airflow and pressure 
ratio at stall as a function of fan blade angle and corrected speed. 
Fan-blade angle is input from the pitch actuator, which is assumed to 
be a simple first order lag with a .1 second time constant. Fan blade 
pitch rate is limited to 100°/ sec. During the RTA simulation study 
(ref . 1) , the pitch actuator was included as part of the power lever 
system model. The actuator model was basically the same as that just 
discussed except that a provision was made for a deadband in blade 
position. During the RTA study the actuator time constant and 
deadband size were varied to investigate the effects of thrust 
response and accuracy setting. 

Compressor corrected airflow and adiabatic efficiency are determined 
as functions of compressor corrected speed (ranging from 65 to 107.5 
percent of design) . 

Fan duct airflow for the turbofan units is computed by subtracting 
compressor airflow from fan airflow. All fan ducts are treated as 
having a fixed percentage total pressure loss. 

Core flow is assumed to be the same as compressor airflow minus 
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compressor bleed. Fuel flow is added at the combustor and the bleed 
airflow is added back into the core flow at the turbines. 

An iteration process is used to compute combustor total pressure PT4 
which is assumed to be equal to compressor discharge static pressure 
PS3. The iteration loop (see ENGNyD subroutines) involves the 
heat-soak lead-lag dynamics and the temperature rise (TT4P-TT3) across 
the combustor due to the fuel flow input WFM. The temperature rise is 
found as a linear function of fuel/air ratio. Total pressure and 
temperature ratios across the high pressure turbine are assumed to be 
constants. Temperature, pressure and flow at the nozzles of the 
turbofan units are calculated assuming mixing of the fan duct and 
low-pressure turbine airstreams. 


Fan, compressor and turbine power are computed from the flow rate 
through the machine times the enthalpy rise or drop across it. 
Enthalpy change is approximated as a constant specific heat times the 
change in temperature. The specific heat values in the equations were 
adjusted to give good agreement with the Detroit Diesel Allison 
steady-state cycle deck representation of the propulsion system. The 
time rate of change of rotor speed is calculated from the difference 
in power absorbed by the compressor or fan and the associated turbine 
power. Power losses due to gearing are neglected. High pressure 
rotor accelerations are calculated individually as shown in figure 4. 
However, since all three fans are connected by shafting, the low 
pressure rotor powers are summed as shown in figure 4. The time rate 
of change of fan speed NLDT is integrated to obtain fan speed which is 
the same for all three fans (see subroutine DVTOL) . An option is 
available to allow declutching of the remote fan. In that case the 
power absorbed by the fan HPF3 and its inertia go to zero through a 
lag term. This feature was not used during simulation of the RTA. 
Another option allows simulation of an engine core failure. In that 
case the power output of the low pressure turbine drops to zero 
instantaneously. In the case of a turbofan failure the temperature of 
the nozzle airflow is set equal to the fan duct air temperature (i.e. 
no heating from core) . 

Thrust calculations are based on conventional momentum equations using 
the appropriate airflows and jet velocities. The nozzle exit 
temperature of the turbofan units is higher than that of the remote 
fan because of mixing with hot air from the core. The higher 
temperature results in higher nozzle velocity and consequently higher 
thrust. In the case of propulsion unit 3, only the fan is assumed to 
produce thrust (core thrust is neglected). 
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Fuel Control Model 


A simplified block diagram of the fuel control ia shown in figure 5. 
The dynamics of the compressor-inlet total-temperature (TT25) sensor 
and the fuel metering valve are implemented as first order lags. The 
time constant of the sensor is a function of airflow W25. Inputs to 
the fuel control are sensed compressor-inlet temperature T25SN 
compressor exit static pressure PS3, high-pressure-rotdr mechanical 
speed NH and demanded high-pressure-rotor corrected speed PCNHRD. 
There are no measurement dynamics associated with PS3 or NH. When 
used with the RTA simulation, the demanded corrected speed was 
determined from the power lever system output. The fuel control 
includes a proportional-plus-integral controller and fuel schedules 
to limit engine acceleration and deceleration. The MIN and MAX and 
limiter blocks shown in the diagram are part of the engine 
accel/decel limit and overtemperature protection. The output of a 
MIN or a MAX block is the smallest or largest of its inputs 
respectively. The control, as programmed, limits mechanical speed 
operation between 57 and 110 percent of design mechanical speed. The 
output of the fuel valve, WFM, goes to the combustor. 


RESULTS AND DISCUSSION 


Steady state results from the propulsion system computer simulation 
were compared to results from a steady-state cycle deck developed by 
Detroit Diesel Allison (DDA) . Both programs were run at fan corrected 
speeds ranging from 70 to 110 percent of design (fan blade pitch angle 
at design) and at fan blade pitch angles ranging from 1 to -12 degrees 
(fan corrected speed of 100%). In the fan corrected speed range of 90 
to 110% the maximum thrust error was 4% and most other variables (e.g. 
fuel flow, turbine inlet temperature) were in error by less than 3%.' 
At 70% fan corrected speed thrust and fuel flow were in error by 7% or 
less and temperatures and pressures were off by less than 3%. 

Figure 6 illustrates the type of transient results that can be 
obtained for internal parameters from the propulsion simulation. Not 
all of the internal engine parameters that are available in the 
simulation are shown. The particular case shown was taken from 
reference 1 and is typical of the results obtained from the RTA 
simulation study. The results are for the last 120 seconds of an 
approach trajectory to hover. Failure of engine number 2 was 
programmed to occur at an altitude of 305m (approximately 30 seconds). 
Cycling of engine parameters such as turbine inlet temperature TT4 is 
evident as is the large jump in TT4 after the engine failure occurs. 
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In fact, at the end of the transient, TT4 approaches the 1 hour 
contingency rating of 1620 K specified in reference 5. Events such 
as these are of interest as to how they impact engine life and design 
criteria for V/STm. a^.raft. Analysis of these data are beyond the 


scope of this report. 


CONCLUDING REMARKS 


A propulsion system model suitable for a non-real-time digital 
simulation of a lift/cruise fan V/STOL aircraft was presented. 
Steady-state agreement with detailed cycle deck calculations is good. 
The model has been integrated with the NASA-Ames mathematical model of 
the V/STOL Research and Technology Aircraft (RTA). It could bo 
adapted _ for use with any similarly configured aircraft. The 
propulsion model is sufficiently detailed to allow investigation of 
thrust response requirements and low-cycle-fatigue/engine-life- 
deterioration during approach and vertical landing trajectories. 
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APPENDIX A 


DEFINITION OF PROPULSION SYSTEM FORTRAN VARIABLES 

The following list defines the basic propulsion system parameters 
including all variables required to interpret the mathematical model 
equations. The general form of the variable names is ABCxy where ABC 
refers to the physical quantity, x refers to the station number (1 or 
2 digits, see fig. 3), and y refers to the propulsion unit number (see 
figs. 1,2). The propulsion unit designation y is attached to the name 
only in subroutine DVTOL. Both SI and English units are given for 
each variable. When used with the RTA simulation, the computations, 
as documented in this report, were made using English units. 

ALTENG aircraft altitude, m (ft) 

BETAF fan blade pitch angle deg 
BETFD or 

BETN command to fan blade actuator 

DELTAT integration step size from LIFAN (typically .01), sec. 

DETETO transient time input from LIFAN (typically .05), sec 
DTT temperature increment above standard atmosphere, K(deg R) 

EIE initial fuel flow rate, Kg/sec (lbm/hr) 

EOKL output of first order lag in LIFAN when fan 3 clutch is 

engaged or disengaged 

ETAC core compressor efficiency, dimensionless 

ETAR inlet pressure recovery, dimensionless 

FGROSS or 

FGROS propulsion unit gross thrust, N (lbf) 

FN propulsion unit net thrust, N(lbf) 

FRAM ram drag, N(lbf) 

FREQ 2 heat soak lag frequency (see fig. 4) rad/sec 

FREQ4 heat soak lag frequency (see fig. 4) rad/sec 

HPF power absorbed by fan, W(hp) 

HP IT power generated by low pressure turbine, W(hp) 

HP2C power absorbed by core compressor, W(hp) 

HP2T power generated by high pressure turbine, W(hp) 

IMODE < 0 for initialization pass 

0 for transient run 
NH compressor speed, rpm 

NHDT compressor acceleration, rpm/sec 

NL fan speed, rpm 

NLDT fan acceleration, rpm/sec 

P » PS static pressure, N/cm**2 (psi) 

PCNHD or 

PCNHRD commanded compressor corrected speed, % of design 
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PCNHR 

PCNLR 

PRSTL 

PS3GS 

PT 

P13Q2 

P25Q2 

SM 

T 

TAUT25 
TAU1 
TAU3 . 
TT • 
TT3P 
TT4P ; 
T13Q2 - 
T25Q2 
T25SN 
V 

VCLUCH 

VJ3 

W 

WA 

WFH 

WFI 

WFM 

WG 

WSTL 

W1R 

W25R 

XFAIL 

XNPCT 


compressor corrected speed, % of design 

fan corrected speed, % of design 

fan stall pressure ratio for given PCNLR 

guess for compressor exit static pressure in iteration loop 
(see ENGNyD subroutines), N/CM**2 (psi) 
total pressure, N/CM**2(psi) 
fan tip total pressure ratio 
fan hub total pressure ratio 
fan stall margin 
; static temperature, K(deg R) 
compressor inlet temperature sensor time constant, sec 
heat soak dynamics time constant (see fig. 4), sec 
heat soak dynamics time constant (see fig. 4), sec 
-total temperature, K (deg R) 

compressor exit total temperature before heat soak K (deg R) 
combustor exit total temperature before heat soak, K (deg R) 
fan tip total temperature ratio 
fan hub total temperature ratio 

output of compressor inlet total temperature sensor, K(deg R) 
flow velocity, m/sec 

signal from LIFAN indicating when fan 3 is engaged (“1 . ) or 
disengaged (=0.) 

propulsion unit 3 low pressure rotor moment of inertia, 
kg-m (slug-ft ) 

inlet and engine core airflow rate, kg/sec (Ibm/sec) 
fan duct airflow rate, kg/sec (lbm/sec) 
fuel flow rate, kg/hr (lbm/hr) 

output of fuel controller integrator, kg/hr (lbm/hr) 
fuel flow rate, kg/sec (lbm/sec) 

mass flow rates downstream of combustor (includes fuel flow), 
kg/sec (lbm/sec) 

fan stall airflow for given PCNLR, kg/sec (lbm/sec) 

fan corrected airflow, kg/sec (lbm/sec) 

compressor corrected airflow, Kg/sec (lbm/sec) 

flag to indicate engine failure (1.0 for engine operating, 0. 

for engine failed) 

fan mechanical speed, %of design 
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APPENDIX B 


DESCRIPTIONS AND LISTINGS OF DIGITAL COMPUTER SUBROUTINES 


Subroutine DVTOL 

This 'is the main controlling routi ne f or, the propulsion system 
computer program. It is assumed that DVTOL is called from a main 
program or another subroutine that supplies the aircraft altitude 
ALTENG andjMach number XMACH as well as initialization constants, 
control flags (e.g. XFAIL) and any other required inputs. References 
to LIFAN in the comment statements of DVTOL refer to input from or 
output to the power lever system part of the RTA airframe model. 
DVTOL sets up the initial conditions and controls the flow of the 
program during the dynamic segment. Output statements have been 
omitted, but just about all fan and engine variables are available 
via the COMMON statements. 
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DTOL ,10/24/80 08^45:47 


100 

200 C 
300 C 
400 
300 
600 
700 
800 
900 
1000 
1100 
1200 
1300 
1400 
1500 
1600 
1700 
1800 
1900 
2000 
2100 
2200 
2300 
2400 
2500 
2600 
2700 
2800 
2900 
3000 
3100 
3200 
3300 
3400 
3500 
3600 
3700 
3800 
3900 
4000 
4100 
4200 
4300 
4400 
4500 
4600 
4700 
4800 
4900 
5000 
5100 
5200 
5300 C 


B 

C 

D 

E 

F 

Q 


SUBROUTINE DVTOL 
AS CONVERTED FOR LIFAN 

COMMON /FMEMR/IX(60),JY(60),IERR(60) 

COMMON /INTVAL/DT1,DT2,IC0UNT 
COMMON /XFLQAT/A(500)/IFIXED/IAC200) 

COMMON /YF10AT/B( 300 ) 

?^nS L i^V0^F5^Si^H2,NH3,NLD T1 ,N LDT 2,H LD T3,NH DT1 ,NH DT 2,NHDT3,MHO 
COMMON /XXI/ ALT,MNO, AE8 , CV8 , CV18 , PS3GS , PT8GS , QX 
COMMON /XX2/XALTC8) ,ZPFA(8) ,2TFA(8),XPCNH2(16),2ETACC16),- 
A ^^MS«tl 1C12 ?/ 2W25R(12),X F 1 8C 1 7) ,Z F1 8(i7) ,XF81(1Q) , 2F81 ( 10 ) , YY( 50 ) ,DYDTC50) 

SSSS /Kb/ ^I T T A ASU;S R R !i:5^^^i ZFrPR "’ 5> ’ 2mR "- 5i 

A c ss5es l r^,s^; z xjii z x5Sk^fi; sHi ' riasNi ' T3sHi * ERLiHi - xi8:i ' T25sH1 '- 

COMMON /DYDT1/ NLDT1,NHDT1 , T3LDT1 , T4LDT1 , PT0DT1, T12DT1 , T3SDT1, ERNDT1 - 
A cSmSSn 1 ;s?DEV- QMVDT 1 ,XMVDT1 ' BETDT1 ' XXADT1 ' XXB ° h ; WFIDTl ;WFH6 t 1 

A A181 , WFH1 > PCNLD1 ,XXKL1 ,XXKP1 > PO 1 ,T01, Pill, Till. TT21»- 

V01,PCNLR1,P13Q21,T13Q21,W11,PT21,PT131,TT131,PT251,TT251,- 
PCNHR1 » W25R1 » W251 ,P13181»PT181» P0Q181 # WA181 ,U31 , W641 » WFM1 >- 
PS31»PT31»P3Q251,TT31,TT41»PT41, TT421 , WG81 , PT81,PT421»“ 

TT81 » HP1T1 » HPF1 , HP2T1 » WFHI1 »HP2C1 »P0Q81 » V81 » V181 » FGROS1 
FRAM1 » FN1 » WFMPS1 , WFMPT1 , WFMTT1 , WFMQK1 , PHI1 , PS3C1 , WFMAC1 , FNZ181, - 
WG181 » WFMDC1 , XMVDC1 * TT181 , ERRNL1 , PCNHD1 , QNH1 , QMV1 , HPX1 » BETFD1 , - 
H SMI 

A c 3R^g 2 ^5$.SIfA^x5!k?x^kfw?II SHa - TiasH2 ' T3SNa ' eRLIM2 ’ ><lsz ’ T25SH2 '- 

A A182,WFH2,PCNLD2,XXKI2,XXKP2,P02,T02,P112,T112,TT22,- 
B V02,PCNLR2,P13Q22,T13Q22,W12,PT22,PT132,TT132,PT252,TT252,- 
p CNHR2,W25R2»W252»P13182,PT182,P0Q182>WA182,W32»WG42>WFM2»- 
PS32»PT32,P3Q252,TT32»TT42»PT42,TT422,WG82,PT82,PT422»- 
TT82»HP1T2»HPF2,HP2T2»WFHI2>HP2C2,P0Q82,V82,V182, FGR0S2, - 
pp AM2, FN 2»WFMPS2,WFMPT2,WFMTT2,WFMQK2,PHI2,PS3C2,WFMAC2,FNZ182,- 
WG182,WFMDC2,XMVDC2,TT182, ERRNL2,PCNHD2,QNH2,QMV2,HPX2,BETFD2,- 
H SM2 

A C QMVLG3^XMV3,BETAF3?XXA3?XXB3?WFI3 SN3 ' T12SN3,T3SN3 ' ERLIN3,X183 ' T25SN3 '~ 
COMMON /SIDE3/- 

A A183 » WFH3 , PCNLD3 ,XXKI3 , XXKP3 , P03,T03,P113,T113»TT23»- 
B V03,PCNLR3»P13Q23,T13Q23,W13,PT23,PT133,TT133,PT253»TT253,- 
PCNHR3,W25R3,W253,P13183,PT183,P0Q183,WA183,W33,WG43,WFM3,- 
PS33,PT33,P3Q253,TT33,TT43,PT43,TT423,WG83,PT83,PT423,- 
TT83,HP1T3»HPF3,HP2T3,WFHI3,HP2C3>P0Q83»V83,V183, FGR0S3 . - 
, F ,^^' FN3 ' UFMPS3 ' WFMPT3 ' WFMTT3 ' WFMQK3 > PH I 3 ' p S3C3,WFMAC3,FNZ183 > - 
WG183,WFMDC3,XMVDC3,TT183,ERRNL3,PCNHD3,QNH3,QMV3,HPX3,BETFD3,- 

H SM3 

ADDED FOR LIFAN CONVERSION 
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DTOL 


,10/24/80 08:45^47 


5400 

5500 

5600 

5700 

5800 

5900 

6000 

6100 

6200 

6300 

6400 

6500 

6600 

6700 

6800 

6900 

7000 

7100 

7200 

7300 

7400 

7500 

7600 

7700 

7800 

7900 

8000 

8100 

8200 

8300 

8400 

8500 

8600 

8700 

8800 

8900 

9000 

9100 

9200 

9300 

9400 

9500 

9600 

9700 

9800 

9900 

10000 

10100 

10200 

10300 

10400 

10500 

10600 


DIMENSION EIE(3),BET(3),DWF(3),DMF(3), FGUC3 ) , 0MEG(3 ) , OMEGFC 3 ) 

DIMENSION PCNHCC3) 

EQUIVALENCE CXMACH , AC7 1 ) ) , C ALTENG, A C83 ) ) , ( D2R, A( 358) ) 

EQUIVALENCE < DELT2 , A( 168 ) ) , ( BET( 1 ) , BC47 ) ) , ( EIE( 1 ) , B (285) > 

EQUIVALENCE (XNPCT, B(121) ) , ( DWF< 1 ) , B( 90 ) ) , (FGUC 1 ) , B( 95) ) 

EQUIVALENCE COMEGC 1 ) , B( 62) ) , ( OMEGFC 1 ) , B( 105 ) ) , ( DETETO, A (168 ) ) 
EQUIVALENCE CDMFCl ), B(87 )), COMEGC 1 ), B( 62) ),( OMEGFC 1 ), BC105) ) 
EQUIVALENCE (XNPCT3, B< 293 )),( EOKL , 5(294 ) ) , CDTPRNT , A( 2) ) , (TIME, AC3 ) ) 
EQUIVALENCE (PCNHC(l) >B(290) ) , (IA( 1) , IMODE) , (PROTIM, A(4) ) , (XFAIL , A(5) ) 
EQUIVALENCE ( BC 120 ) , VCLUCH ) , ( BC294 > , EOKL ) , (DELTAT , AC 16 9) > , (DTT , A (170 ) ) 
EQUIVALENCE (ETAR, A(6>) 

C 

IF(IMODE.GE.O) GO TO 233 
DO 2 1=1,60 
IX(I)=1 
JY( I) =1 
IERRC I ) =1 
2 CONTINUE 
CKKMKKINITIAL CONDITIONS 
C INPUT FOR FAN BLADE PITCH FROM LIFAN 
YY( 1 ) =3600 • 

YY<2)=14528. 

YY( 3 ) =1221 • 

YY(4)=2493. 

YYC5)=EIEC1) 

YYC6 )=EIE(2) 

YY( 7 ) =EIE( 3 ) 

YY( 8 ) =0 . 0 
YY( 9 ) -0 . 0 
YY< 10 ) =583 . 8 
YY(11)=0.0 
YYC 12) =0 . 0 
YY( 13 ) =BETAF1 
YY( 14 ) =0 . Q 
YY<15)=3600. 

YY( 16 ) =14528 . 

YY ( 17 ) = 1221 . 

YY( 18 ) =2493 . 

C FUEL FLOW INPUT FROM LAST WFH 
YYC 19) =EIE( 1 ) 

YYC20 ) =EIE(2 ) 

YYC21 )=EIEC3) 

C 

C 

YY(22)=0.0 
YYC 23 ) =0 . 0 
YYC24 ) =583 . 8 
YY( 25 ) = 0 . 0 
YY(26)=0.0 
YYC27 ) =BETAF2 
YYC28) -0 , 0 
YYC29)=0 .0 
YY( 30 ) =3600 . 
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DTOL 


,10/24/80 08 * 45 * 47 


10700 

10800 

10900 

11000 

11100 

11200 

11300 

11400 

11500 

11600 

11700 

11800 

11900 

12000 

12100 

12200 

12300 

12400 

12500 

12600 

12700 

12800 

12900 

13000 

13100 

13200 

13300 

13400 

13500 

13600 

13700 

13800 

13900 

14000 

14100 

14200 

14300 

14400 

14500 

14600 

14700 

14800 

14900 

15000 

15100 

15200 

15300 

15400 

15500 

15600 

15700 

15800 

15900 


YY(31)=14576. 

YY(32)=1197. 

YY<33)=2539. 

YY(34)=0 • 0 

YY<35)=0.0 

YY( 36 ) =0 . 0 

YY(37)=0.0 

YYC38)=0.0 

YY<39)=549.7 

YY(40)=0.0 

YY(41 ) =0 . 0 

YY(42)=BETAF3 

YY(43)=0 . 0 

YY(44)=0 . 0 

CxxxxxTIMING CONSTANTS 
DT1=DELTAT 
DT2=0.5XDELTAT 
TIME=0.0 
233 CONTINUE 
UFHSV1=0. 

WFHSV2=0. 

WFHSV3=0. 

KFUEL=0 

CXXXXXMACH NUMBER INPUT FROM LIFAN 
C 

MNO=XMACH 

CXXXXXALTITUDE INPUT FROM LIFAN 
C 

ALT=ALTENG 

CXXXXXENGINE ENVIRONMENTAL CALCULATIONS 
P01=FUN1C1,8,XALT,ZPFA,ALT) 
T01=FUN1(2,8,XALT,ZTFA,ALT) 
CXXXXXNONSTANDARD DAY TEMP. 

T01=T01+DTT 
T02=T01 
T03=T0 1 
P02=P01 
P03=P01 

CXXXXXINLET CALCULATIONS 

PT21=ETARXP01X(1.+0.2XMNOXX2)XX3.5 

TT21=T01X(1.+0.2XMNOXX2) 

PT22=PT21 

PT23=PT21 

TT22=TT21 

TT23=TT21 

CXXKXXBET AND PCNHC INPUT FROM LIFAN 
BETAF1=BET ( 1 ) 

BETAF2=BET (2 ) 

BETAF3=BET(3) 

PCNHDl=13.89X(PCNHC(l))XX.2+28.87 
PCNHD2=1 3.8 9X(PCNHC( 2) )xx. 2+28.87 
PCNHD3=17.07X(PCNHC(3))XX.2+16.11 
VJ3=VCLUCHX14.4 
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DtOL 
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16000 

16100 

16200 

16300 

16400 

16500 

16600 

16700 

16800 

16900 

17000 

17100 

17200 

17300 

17400 

17500 

17600 

17700 

17800 

17900 

18000 

18100 

18200 

18300 

18400 

18500 

18600 

18700 

18800 

18900 

19000 

19100 

19200 

19300 

19400 

19500 

19600 

19700 

19800 

19900 

20000 

20100 

20200 

20300 

20400 

20500 

20600 

20700 

20800 

20900 

21000 

21100 

21200 


BETFD1=BETAF1 

BETFD2=BETAF2 

BETFD3=BETAF3 

C SET RUNTIM a TO DT2 FROM LIFAN 
RUNTIM = DETETO 
IQUIT=RUNTIM/DELTAT+.l 

C WRITEC 6 >101) 

CXXKKXDYNAMIC SEGMENT OF MODEL 
DO 700 LOOPl=l,IQUIT 
ICOUNT=-l 

600 CONTINUE 
NL1=YYC1) 

NH1=YY(2) 

T3L1=YYC3> 

T4L1=YYC4) 

WFI1=YY(5) 

WFI2=YY(6) 

WFI3 S YY(7 ) 

ERLIN1=YYC8) 

TWIST=YY(9) 

T25SN1=YYC10) 

QMVLG1=YY(11) 

XMV1=YY(12) 

DUM14=YY(14) 

NL2=YY( 15 ) 

NH2=YY( 16 ) 

T3L2=YY(17) 

T4L2=YY(18) 

WFH1=YY(19) 

WFH2=YYC20) 

WFH3=YY(21) 

ERLIN2=YY(22) 

DUM24=YYC23) 

T25SN2=YY(24) 

QMVLG2=YY(25) 

XMV2=YYC26) 

DUM25=YY( 28) 

DUM26=YY(29) 

NL3=YY( 30 ) 

NH3=YY( 31 ) 

T3L3=YY(32) 

T4L3=YY(33) 

DUM31=YY<34) 

DUM32=YY(35) 

DUM33=YY ( 36 ) 

ERL IN3=YY ( 37 ) 

DUM34=YYC38) 

T25SN3=YY(39) 

QMVLG3=YY (40 ) 

XMV3=YYC41) 

DUM35=YY ( 43 ) 

DUM36=YY(44) 

CALL ENGN1D 
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DTOL 
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21300 

21400 

21500 

21600 

21700 

21800 . 

21900 

22000 

22100 

22200 

22300 

22400 

22500 

22600 

22700 

22800 

22900 

23000 

23100 

23200 

23300 

23400 

23500 

23600 

23700 

23800 

23900 

24000 

24100 

24200 

24300 

24400 

24500 

24600 

24700 

24800 

24900 

25000 

25100 

25200 

25300 

25400 

25500 

25600 

25700 

25800 

25900 

26000 

26100 

26200 

26300 

26400 

26500 


CALL ENGN2D 
CALL ENGN3D 

IFdMODE.LT, 0) GO TO 216 
ERNDTl=PCNLDl-NLl/36. 

ERNDT2=ERNDT1 

ERNDT3=ERNDT1 

HP1T1=HP1T1*XFAIL 

HP1TT=HP1T1+HP1T2+HP1T3 

HPFT=HPF1+HPF2+HPF3 

DELHP=HP1TT-HPFT 

NLDT1=(50154./(2.K14.4+VJ3))*DELHP/NL1 

IFC ABSCNLDT1 ) . GT , 9549.3) NLDT1 = SIGN(2549 . 3> NLDT1 ) 

NLDT2=NLDT1 

NLDT3=NLDT1 

DYDT(1)=NLDT1 

DYDT(2)=NHDT1 

DYDT(3)=T3LDT1 

DYDT (4)=T4LDT1 

DYDT ( 5 ) =WFIDT1 

DYDT ( 6 ) =WFIDT2 

DYDT ( 7 ) = WFIDT3 

DYDTC8)=ERHDT1 

DYDT C 9 ) =0 . 0 

DYDT(10)=T25DT1 

DYDT(11)=QMVDT1 

DYDT C 12 ) =XMVDT1 

DYDT ( 13 ) =0 . 

DYDT ( 14 ) =0 . 0 
DYDT ( 15) =NLDT2 
DYDTC16 )=NHDT2 
DYDTC17)=T3LDT2 
DYDT ( 18 ) =T4LDT2 
DYDT(19)=WFHDT1 
DYDT(2Q)=WFHDT2 
DYDT (21 ) =WFHDT3 
DYDT (22)=ERNDT2 
DYDT (23 ) =0 . 0 
DYDT(24)=T25DT2 
DYDT(25)=QMVDT2 
DYDT ( 26 ) =XMVDT2 
DYDT (27 ) = 0 . 

DYDT (28 ) =0 . 0 
DYDT ( 29 ) =0 . 0 
DYDT ( 30 ) =NLDT3 
DYDT ( 31 ) =NHDT3 
DYDTC32)=T3LDT3 
DYDT (33 ) =T4LDT3 
DYDT ( 34 ) =0 . 0 
DYDT ( 35 ) =0 . 0 
DYDT ( 36 ) =0 , 0 
DYDTC37 ) =ERNDT3 
DYDT ( 38 ) =0 . 0 
DYDT ( 39 ) =T25DT3 


15 



DTOL 
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26600 
26700 
26800 
26900 
27000 
27100 
27200 
27300 
27400 
27500 
27600 
27700 
27800 
27900 
28000 
28100 
28200 
28300 
28400 
28500 
28600 
28700 
28800 
28900 
29000 
29100 
29200 
29300 
29400 
29500 
29600 
29700 
29800 
29900 
30000 
30100 
30200 
30300 
30400 
30500 
306 00 
30700 
30800 
30900 
31000 
31100 


DYDTC40 )=QMVDT3 
DYDT(41)=XMVDT3 
DYDT (42 ) =0 . 

DYDT ( 43 ) =0 . 0 
DYDT C 44 ) =0 . 0 

Cxxxxh INTEGRATE DIFFERENTIALS 
CALL EULERM(DYDT, YY) 
YY(30)=YY(30)*EOKL 
IF CICOUNT) 600,610,650 
650 CONTINUE 

KFUEL=KFUEL+1 
WFHSV1=WFH1+WFHSV1 
WFHSV2=WFH2+WFHSV2 
WFHSV3=WFH3+WFHSV3 
CxxxxxUPDATE TIME 

TIME=TIME+DELTAT 
GO TO 700 
610 CONTINUE 
GO TO 600 
700 CONTINUE 

C OUTPUT VALUES FOR LIFAN 
WFHOU1 = WFHSV1/KFUEL 
WFH0U2 = WFHSV2/KFUEL 
WFH0U3 = WFHSV3/KFUEL 
216 XNPCT=NLl/36. 
XNPCT3=NL3/36. 

DWF( 1 ) =WFH0U1 
DWF(2)=WFH0U2 
DWF(3)=WFH0U3 
FGUC1)=FGR0S1 
FGU(2)=FGR0S2 
FGU ( 3 ) =FGR0S3 
DMF(1)=WA181X. 0310559 
DMF(2)=WA182x .0310559 
DMF(3)=WA183X. 0310559 
0MEG(1)=NH1X. 1047197 
0MEG(2)=NH2X. 1047197 
0MEG(3)=NH3*. 1047197 
0MEGF(1)=NL1X. 1047197 
0MEGF(2)=NL2X. 1047197 
OMEGFC3 ) =NL3X .1047197 
EIE(1)=WFH1 
EIEC2)=WFH2 
EIE(3)=WFH3 
215 RETURN 
END 
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Subroutine ENGN1D 


This subroutine Is called from DVTOL and computes variables for 
turbofan propulsion unit number one (see fig, 1). Included are fan, 
core, core inlet temperature Bensor, core speed (fuel) control, and 
fan-blade pitch actuator parameters. (It should be recallod that the 
pitch actuator waB handled outside the propulsion model for the RTA 
simulation.) Derivative terms are also calculated except for fan 
acceleration NLDT, which requires inputs from all three propulsion 
units. 
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100 

200 

300 

400 

500 

600 

700 

800 

900 

1000 

1100 

1200 

1300 

1400 

1500 

1600 

1700 

1800 

1900 

2000 

2100 

2200 

2300 

2400 

2500 

2600 

2700 

2800 

2900 

3000 

3100 

3200 

3300 

3400 

3500 

3600 

3700 

3800 

3900 

4000 

4100 

4200 

4300 

4400 

4500 

4600 

4700 

4800 

4900 

5000 

5100 

5200 

5300 


SUBROUTINE ENGN1D 

REAL MNO,NL,NLDT,NH,NHDT,NHD,NHERR 
COMMON /XXI/ ALT,MN0, AE8 , CV8 , CV18 , PS3GS , PT8GS , QX 
COMMON /XX2/ XALT(8),ZPFA(8),ZTFA(8),XPCNH2(16),ZETACC16),- 
A XPCNH1(12),ZW25R(12),XF18(17),ZF18(17),XF81(10),ZF81(10),YY(50),DYDTC50) 
COMMON /FAN/ XBETA ( 7 ) , YPCNLR ( 5 ) , ZFFLOWC 7 , 5 ) ,ZFTPR ( 7 , 5 ) , ZFTTRC 7 , 5) 

COMMON /HUB/ XBETA1 (8 ) , ZFHPR (8 , 5 ) , ZFHTRC8 , 5 ) 

C0MM0N/C0RE/XNH(24),YTT25(5) ,ZWFSCH(24 , 5 ) 
C0MM0N/FAN3TL/XBETA2(5),ZWSTL(5,5),ZPRSTl(5,5) 

COMMON /YY1/ NL,NH,T3L,T4L,PT0SN,T12SN,T3SN,ERLINT,X18,T25SN,QMVLG,XMV,- 
A BETAF, XXA>XXB ,WFI 
COMMON /XFLOAT/AC500) 

EQUIVALENCE CXFAIL,A(5>) 

COMMON /DYDT1/ NLDT , NHDT , T3LDT , T4LDT , PTOSDT, T12SDT , T3SDT, ERNLDT,X18DT, - 
A TT25DT ,QMVLDT ,XMVDT ,BETADT »XXADT ,XXBDT,WFIDT, WFHDT 
COMMON /SIDE1/- 


A 

B 

C 

D 

E 

F 

G 


A18 

VO 


, WFH 

, PCNLRD,XXKI , XXKP 

,P0 

, TO 

,P11 

, Til 

, TT2 

, PCNLR 

, P13Q2 , T13Q2 ,W1 

, PT2 

, PT13 

, TT13 

, PT25 

t TT25 

, LI25R 

, W25 ,P13Q18,PT18 

, P0Q18 

» WA18 

,W3 

, WG4 

,WFM 

, PT3 

>P3QP25,TT3 ,TT4 

, PT4 

,TT42 

, WG8 

,PT8 

, PT42 

, HP1T 

, HPF , HP2T , WFHI 

, HP2C 

, P0Q8 

,V 8 

, V18 

,FGROSS,- 

, FN 

,WFMPSS,WFMPTP , WFMPST 

, WFMQK 

, PHI 

,PS3C 

, WFACC, 

FNZl8 

, WFDEC, 

XMVDEC, TT18 , ERRNL , 

PCNHRD, 

QNH 

, QMV 

, HPX 

BETAFD, - 


PS3 
TT8 
FRAM 
WG18 
H SM 
ITR3S=0 

RTH2=SQRT( TT2/518 . 67 ) 
QRTH2=1./RTH2 
RTT0=SQRT (TO ) 

V0=49 . 018*MN0xRTT0 
CxxxxxENGINE MODEL 
600 CONTINUE 
CxxxxxFAN CALCULATIONS 

PCNLR=NL/36./(SQRT(TT 2/556. ) ) 


W1R=FUN2 ( 10, 7, 5, XBETA, YPCNLR, ZFFLOW, BETAF, PCNLR) 
P13Q2=FUN2( I 1,7, 5, XBETA, YPCNLR ,ZFTPR, BETAF, PCNLR) 
T13Q2=FUN2 (12, 7, 5, XBETA, YPCNLR, ZFTTR, BETAF, PCNLR) 

IFCWIR .LT. 0.) W1R=0. 

IF(T13Q2 .LT. 1.) T13Q2=1. 

IF(P13Q2 .LT. 1.) P13Q2=1. 

W1=W1RX(PT2/1A.696)XQRTH2 
CHKHXKFAN STALL MARGIN CALCULATIONS 

WSTL=FUN2(4 , 5 , 5 » XBETA2, YPCNLR, ZWSTL , BETAF, PCNLR) 

PRSTL=FUN2(5,5,5,XBETA2,YPCNLR,ZPRSTL, BETAF, PCNLR) 

SM=100 ,x( (PRSTLXWlR/pi3Q2/(WSTLX . 0 929X17 . 18/. 4536 ) )-l . ) 

P25Q2=FUN2(13,8,5,XBETA1, YPCNLR, ZFHPR, BETAF, PCNLR) 

T25Q2=FUN2 (14, 8, 5, XBETA1, YPCNLR, ZFHTR, BETAF, PCNLR) 

PT13=P13Q2XPT2 

TT13=T13Q2XTT2 

PT25=P25Q2XPT2 

TT25=T25Q2XTT2 


QRTH25=1.0/SQRT(TT25X1.928E-3) 
CXXXXXCOMPRESSOR CALCULATIONS 

PCNHR=NH/154.5/(SQRT(TT25/605.4)) 
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5400 
5500 
5600 
5700 
5800 
5900 
6000 
6100 
6200 
6300 
6400 
6500 
6600 
6700 
6800 
6900 
7000 
7100 
7200 
7300 
7400 
7500 
7600 
7700 
7800 
7900 
8000 
8100 
8200 
8300 
8400 
8500 
8600 
8700 
8800 
8900 
9000 
9100 
9200 
9300 
9400 
9500 
9600 
9700 
9800 
9900 
1 0000 
10100 
10200 
1 0300 
10400 
10500 
1 06 00 


c?Ar:^S}^r U ^ 1C15 ' 12 ^ XPCNH1 ' 2W25R >PCNHR) 

ETAC-FUN1C16, 16 , XPCNH2 , ZETAC , PCNHR) 

W25=W25RX. 068046KPT25XQRTH25 
CxxxxxFAN DUCT CALCULATIONS 
PT18= . 978XPT13 
WA18=W1-W25 
34 W3= . 9036XW25 
WFM=WFH/3600. 

WG4=WFM+W3+. 025XW25 

FREQ2=0.0104XW3 

TAU1=0.843/FREQ2 

FREQ4=0.0240XWG4 

TAU3=0.868/FREQ4 

s?!! or disch#ri!e sutic pressure 

KPS3 = 0 

33 ITR3S=ITR3S+1 
PS3=PS3GS 

IFCPS3.LT. 1.0) PS3=1 . 0 
PT3=1 . 057XPS3 
P3QP25=PT3/PT25 

T3Q25=1.+(P3QP25XX. 2857-1. )/ETAC 

TT3P=T3Q25XTT25 

T3LDT=(TT3P-T3L)XFREQ2 

TT3=TAU1XT3LDT+T3L 

TT4P=53233.#(WFM/W3)+159.+TT3 

T4LDT=CTT4P-T4L)XFREQ4 

TT4=TAU3XT4LDT+T4|. 

PS3=. 0812XWQ4XSQRTCTT4) 

ERR=PS3-PS3GS 

IFCKPSE^EqIi )G 0 R T0 P f 0 

GO TO 35 
40 CONTINUE 
PT4=PS3 

CxxxxxHIGH AND LOW TURBINE EXIT CONDITIONS 
TT42= . 7897XTT4 
WG8=W25+WFM 
PT42=.300XPT4 
PT8=. 01557MWG8+. 919XPT18 
TT8Q42=1.-(1. -CPT8/PT42) XX . 248 )X .866 
TT8=TT42XTT8Q42 
IFCXFAIL .LT. 0.1) TT8=TT25 
CxxxxxFAN ROTOR DYNAMICS 

HP1T=.3817XWG8X(TT42-TT8) 

HPF= . 3356XCW25X(TT25-TT2)+WA18X(TT13-TT2) ) 

CXXXXXCORE ROTOR DYNAMICS 

HP2T=.417XWG4X(TT4-TT42) 

HP2C=.369XW3x(TT3-TT25) 

NHDT =37710. /NHXCHP2T-HP2C) 

CXXXXXTHRUST CALCULATIONS 
WG18=WG8+WA18 

TT18=(WA18XTT13+WG8XTT8)/WG18 
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10700 

10800 

10900 

11000 

11100 

11200 

11300 

11400 

11500 

11600 

11700 

11800 

11900 

12000 

12100 

12200 

12300 

12400 

12500 

12600 

12700 

12800 

12900 

13000 

13100 

13200 

13300 

13400 

13500 

13600 

13700 

13800 

13900 

14000 

14100 

14200 

14300 

14400 

14500 

14600 

14700 

14800 


P0Q18=PG/PT18 

FNZ18=FUN1(17,17,XF18,ZF18,P0Q18) 

CxhkhhoR SUBSTITUTE FNZ18=SQRT( 1 . -PQQ18** (2 ./7 . ) ) 

V18=109.6*SQRT(TT18)*FNZ18XCV18 

FGROSS=WG18XV18*.0311 
54 FRAM=W1*V0K0 . 031085 
FN=FGRGSS-FRAM 

CXXXXXCORE INLET AIR TEMPERATURE SENSOR 

TAUT25=FUN1(18,10,XF81,ZF81,1.38*W25) 

TT25DT=(TT25-T25SN)/TAUT25 

TH25SN=T25SN/605.4 

RTH25S=SQRTCTH25SN) 

CXXXXXFUEL CONTROL 

NHD=154.5*PCNHRD*RTH25S 
NHD=AMIN1(NHD, 17000. ) 

NHD=AMAX1(NHD,8795. ) 

NHERR=NH-NHD 

CONST=SQRT(UFH)HNH*1.494E-6 
WFPR0P=-2 .29*NHERR*C0NST 
CON=-3 . 21 

IFCWFC2 . NE. WFGV) CON=0. 

WFIDT=NHERR*CON*CONST 

WFGOV=WFPROP+WFI 

WFSCH=FUN2(19,24,5,XNH,YTT25,ZWFSCH,NH,T25SN) 

WFACMX=- .0134XNH+234.1 

IFCWFSCH.GT. WFACMX) WFSCH=WFACMX 

IFCPS3 • LE. 200 . ) PC=PS3 

I F( PS3 . GT • 200 • ) PC=200.+(PS3-200. )*C-2. ) 

WFACC=WFSCH*PC 

WFDEC = , 385XWFACC 

WFC1=AMIN1(WFACC»WFG0V) 

WFC1=AMAX1(WFC1,WFDEC) 

WFC1=AMIN1(WFC1,4438. ) 

WFC1 = AMAX1 CWFC1 #295*858) 

WFC2=WFC1 

WFGV=WFGOV 

WFHDT=(WFCl-WFH)/.05 
CXXXXXFAN PITCH ACTUATOR 

BETADT=10 .X(BETAFD-BETAF) 

IFCABS(BETADT) .GT.1G0. ) BETADT=SIGNC 100 . , BETADT) 

RETURN 

END 
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Subroutine ENGN2D 


This subroutine is called from DVTOL and computes variables for 
turbofan propulsion unit number 2. It is essentially the same as 
ENGN1D except that it was not programmed to include the core failure 
option XFAIL of ENGN1D. This feature could easily be added. 
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SUBROUTINE ENGN20 

REAL MNO , NL > NLDT, NH, NHDT , NHD, NHERR 

COMMON /XXI/ ALT,MN0,AE8,CV8,CV18,PS3GS»PT8GS,QX 

COMMON /XX2/ XALT(8 ) »ZPFA(8) I ZTFA(8),XPCNH2( 16) 7FTAr*M A 1 

C0MM0N/C0RE/XNHC24 ) , YTT25 ( 5 ) ZWFSCHC24 ? 5) 

COMMON/ FANSTL/XB ET A2(5 > , ZWSTL (5 ,5) VZPRSTl(5,5) 

A C BETAF,XXA?XXBtwF^ T3L ' T4L ' PT0SN ' T12SN ' T3SH ' ERLINT ' X18 ' T25SN ' QMVLG ' XMV '” 

,P0 
,PT2 


A18 

VO 


, WFH 
, PCNLR 


PCNHR > W25R 
, PT3 
* HP1T 
» FN 


,T0 
» PT13 


, P0Q18 , WA18 


,PCNLRD,XXKI , XXKP 
, P13Q2 , T13Q2 ,W1 
,U25 , P13Q18 , PT18 

» P3QP25 , TT3 ,TT4 
* HPF » HP2T , WFHI ,... „„ ,, u „ 

,WFMPSS,WFMPTP,WFMPST,WFMQK , PHI 
> WFDEC»XMVDEC» TT18 , ERRNL ,PCNHRD,QNH 


,PT4 
, HP2C 


, TT42 
,P0Q8 


,P11 
# TT13 
,W3 
»WG8 
,V 8 
» PS3C 
, QMV 


PS3 
TT8 
FRAM 
WG18 
H SM 
ITR3S=0 

RTH2=SQRTCTT2/518.67) 

QRTH2=1./RTH2 

RTT0=SQRT(T0) 

VO =4 9. 018XMN0XRTT0 
CXXXXXENGINE MODEL 
600 CONTINUE 
CxxxxxFAN CALCULATIONS 

PCNLR=NL/36./(SQRT(TT2/556. )) 

P^fi5-cMMi?i?'?'c B ^I^;T P S NLR ' 2FFL0M ' BETAF ' PCN >-R) 

P13Q2-FUN2(21»7,5iXBETA,YPCNLR,ZFTPR,BETAF,PCNLR) 
IF(Wir U Lf <2 o')'wiR = 0 TA,YPCNLR,ZFTTR>DETAF ' PC ^^ 

IF(T13Q2 . LT. * 1 . ) T13Q2=1. 

IF(P13Q2 .LT. 1.) P13Q2=1. 

W1=W1RX(PT2/14.696)XQRTH2 
CXXXXXFAN STALL MARGION CALCULATIONS 

WSTL=FUN2(6,5,5,XBETA2,YPCNLR,ZWSTL,BETAF, PCNLR) 

PRSTL = FUN2(7,5,5,XBETA2,YPCNLR,ZPRSTUBETAF,PCNLR) 

SM=100.XC(PRSTLXW1R/P13Q 2 /(WSTLX. 0929*17. 18/ . 4536 ) )-l . ) 
P25Q2-FUN2(23,8»5,XBETA1> YPCNLR,ZFHPR, BETAF, PCNLR) 

PT13=P33(>2PPT2 A ^ ^ 

TT13=T13Q2XTT2 

PT25=P25Q2XPT2 

TT25=T25Q2XTT2 

QRTH25=1.0/SQRT(TT25*1.928E-3) 

CxxkxxcOMPRESSOR CALCULATIONS 

PCNHR=NH/154.5/(SQRT(TT25/605.4)) 

c?^ic.mwor U t , } ( S 5 ' 12,XPCNH1>ZW25R ' PCN HR) 

ETAC-FUN1C26 » 16 ,XPCNH2,ZETAC, PCNHR) 


> Til , TT2 
» PT25 , TT25 
» WG4 ,WFM 
, PT8 > PT42 
» V18 , FGROSS 

, WFACC, FNZ18 
» HPX ,BETAFD,- 
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W25=W25R#.068046*PT25KQRTH25 
C**K#*FAN DUCT CALCULATIONS 
PT18=.978*PT13 
WA18=W1-W25 

34 W3=.9036#W25 
WFM=WFH/36 00 . 

WG4=WFM+W3+.025#W25 
FREQ2=0.0104*W3 
TAU1=Q.843/FREQ2 
FREQ4=0 . 0240KWG4 
TAU3=0.868/FREQ4 

CKHKHKITERATE FOR COMPRESSOR DISCHARGE STATIC PRESSURE 
CKKKHKKKHEAT SOAK DYNAMICS 
KPS3=0 

35 ITR3S=ITR3S+1 
PS3=PS3GS 

IFCPS3.LT. 1.0) PS3=1 . 0 
PT3 = 1 . 057XP53 
P3QP25=PT3/PT25 

T3Q25=1. +CP3QP25**. 2857-1. )/ETAC 

TT3P=T3Q25HTT25 

T3LDT=(TT3P-T3L)*FREQ2 

TT3=TAU1*T3LDT+T3L 

TT4P =53233. X(WFM/W3)+ 159. +TT3 

T4LDT=CTT4P-T4L)*FREQ4 

TT4=TAU3*T4LDT+T4L 

PS3=.0812*WG4*SQRT(TT4) 

ERR=PS3-PS3GS 

CALL CNVRG1C5,ERR,PS3,PS3G$, . 0002, ITR3S,KPS3) 
IFCKPS3 . EQ . 1 )GO TO 40 
GO TO 35 
40 CONTINUE 
PT4=PS3 

CXXXXXHIGH AND LOW TURBINE EXIT CONDITIONS 
TT42= . 7897XTT4 
WG8=W25+WFM 
PT42=.300*PT4 
PT8=. 01557 *WG8+. 91 9*PT18 
TT8Q42=1.-C1 .-(PT8/PT42)**,248)X.866 
TT8=TT42*TT8Q42 
CXXXXXFAN ROTOR DYNAMICS 

HP1T=.3817*WG8K(TT42-TT8) 
HPF=.3356X(W25#CTT25-TT2)+WA18*(TT13-TT2)) 
CXXXXXCORE ROTOR DYNAMICS 

HP2T=.417XWG4*(TT4-TT42) 

HP2C=.369*W3*(TT3-TT25) 

NHDT=37710 ./NH*CHP2T-HP2C) 

C*xxx*THRUST CALCULATIONS 
WG18=WG8+WA18 

TT18=(WA18*TT13+WG8*TT8)/WG18 

P0Q18=P0/PT18 

FNZ18=FUN1(27,17,XF18,ZF18,P0Q18) 

CXXXXXOR SUBSTITUTE FNZ18=SQRT( 1 . -P0Q18**C2 ./7 . ) ) 
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10700 
10800 
1 C 9 0 0 
11000 
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] 18 0C 
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; 2 l 0 ft 
3 2200 
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t C r « r 
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•2/: GO 
12900 
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1 2 1 0 i; 
1 2 2 0 0 
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5 36 C Q 
13700 

1 j o 0 0 
} 3 7: 0 " 

2 9 0 0 0 

1 <t 1 0 0 


:.£*90 


' l'. 1 * 


V18=109-6KSQRTCTT15)«FH216«CV1S 
FGROSS=W018MV18».03U 
5*+ FRAM=W1>?V0H0. 031085 
FM=FGR05S“FRAM 

CHHHKHCORE INLET AIR TEMPERATURE SENSOR 
T AUT25 = FUN1 ( 28 , 1 0 , XF8 1 , ZFS 1 f 1 . 

TT25DT=( TT25-T25SH )/TAUT23 

TH25SN=T25SN/6G5.A 

RTH25S=SQRT(TH25SH) 

C«>Hf^^FUEL CONTROL 

NHD=154 . 5«PCNHRD*RTH25S 
MHD-AMINKHHD/ 17000 . ) 

HHD=AMAX1(NHD,8795. ) 

MHERR=NH-HHD 

C0NST=SQRT(WFH)«KH5<1 .4 94E--0 
W r P R 0 F = - 2 . 2 9 K N H E R R * C 0 M S T 
CON=~3 . 21 

IFCWFC2 . HE.MFGV) COH=Q. 

WFIDT=NKERRKCOHKCONST 

LiFGOV=WFPROP+WFI 

LJFSCH = FUN2(29,24,5,XNH,YT I 23 , Zl !F5CiL XI T*25 3iJ) 

WFACMX=- . 013CKHH+234 . 1 

IFCL4FSCH . GT . WFACMX) WFSCH=!-JFACMX 

IFCPG3 . LE.200 . ) PC-PS3 

IFCPS3.GT.2Q0. ) PC=200 . + (FS3-20Q 

HFACC=WFSCHHPC 

WFDEC=.335*UFACC 

1 1 F C 1 = A M I N 1 ( L ’ F A C C , H F G 0 V ) 

WFC1=AMAX1(UFC1,WFDEC) 

W F C 1 = A M I N 1 C M F C 1 , 4 4 3 8 . > 

W FC 1 = AM AXK W FC 1 , 2 9 5 . 8 58 > 

WFC2=WFC1 

WFGV=WFGOV 

WFHDT = CWFC1-WFH)/ . 05 
CiWMXFAN PITCH ACTUATOR 

E ET A DT = 1 0 . x ( 3 ET A FD - b E 7 A F ) 

I FC ABS (BET ADT > , OT . 1G 0 . ) l z \ /,r?T - 2 3 OX < 1 0 1 . « El APT > 

RETURN 

END 


Subroutine ENGN3D 


This subroutine is called from DVTOL and computes variables for the 
turbine engine and remote fan propulsion unit number 3. The equations 
of ENGN1D have been modified to allow separation of the fan and 
turbojet units • ENGN3D does not contain the core failure option but 
could be easily modified to do so. 
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SUBROUTINE ENGN3D 

REAL MNO,NL,NLDT,NH,NHDT,NHD,NHERR 
COMMON /XXI/ ALT»MN0 » AE8, CV8 , CV18 , PS3GS, PT8GS , QX 
COMMON /XX2/ XALT(8),ZPFA(8),ZTFA(8),XPCNH2(16) , ZETACC 16 ) » - 
A XPCNH1 (12) ,ZI425R(12) ,XF18(17),ZF18(17 ) ,XF81 (10),ZF81(10),YY(50), DYDTC50 ) 
COMMON /FAN/ XBETAC7 ) , YPCNLRC5) ,ZFFL0WC7 , 5) ,ZFTPR(7,5 ,ZFTTR (7,5) 
C0MM0N/C0RE/XNH(24) »YTT25(5) »Z14FSCH(24,5) 

C0MM0N/FANSTL/XBETA2C5) ,ZWSTL(5,5) ,ZPRSTL(5,5) 

A C BETAF > XXA!xXB^MFl' T3UT4UPT ° SN ' T12SN ’ T3SN,ERLINT ' X18 ' T25SN ' QMVLG ' Xr ' 1V '" 

COMMOn'/DYDT3/'nLDT,NHDT,T3LDT,T4LDT,PTOSDT,T12SDT,T3SDT,ERNLDT,X18DT,- 

A *IJ222 T,QMVLDT ' XMVDT ' BETADT,XXADT ' XXBDT ' WFIDT ' WF HDT 

COMMON /SIDE3/- 

,P0 
> PT2 
> P0Q18 
>PT4 
» HP2C 


A A18 
B VO 


PCNHR ,1425R 


»14FH >PCNLRD,XXKI 
*PCNLR » P13Q2 , T13Q2 


>XXKP 
,141 


,1425 ,P13Q18*PT18 

,P3QP25,TT3 ,TT4 
»HPF , HP2T , WFHI 


, TO 
, PT13 


, PT3 

>HP1T 

,FN ,WFMPSS,I4FMPTP,14FMPST,14FM<)K ,PHI 
,WFDEC,XMVDEC, TT18 , ERRNL ,PCNHRD,QNH 


»WA18 
> TT42 
, P0Q8 


,P11 
, TT13 
,143 
,I4G8 
,V8 
, PS3C 
, QMV 


PS3 
TT8 
FRAM 
14G18 
H SM 
ITR3S S 0 

RTH2=SQRT(TT2/318.67) 

QRTH2=1./RTH2 

RTTOsSQRT(TO) 

V0=49.018*MN0XRTT0 
CkxmxkeNGINE MODEL 
600 CONTINUE 
C*xxx#FAN CALCULATIONS 

PCNLR=NL/36 ./CSQRT(TT2/556 . ) ) 

IF(PCNLR .GT. 0.) GO TO 1 
W1R=0. 

P13Q2=1. 

T13Q2=1. 

GO TO 2 

1 W1R=FUN2(30 , 7,5,XBETA, YPCNLR,ZFFL 014, BETA F, PCNLR) 
P13Q2 = FUN2(31,7,5,XBETA> YPCNLR,ZFTPR, BET AF, PCNLR) 

Ti?. Q ,?« FU J l S C 32 'T' 5,XBETA ' YPCNLR » ZFTTR ' BETAF, PCNLR) 

IF(W1R .LT. 0.) W1R=0. 

IF(T13Q2 .LT. 1.) T13Q2=1.0 
IF(P13Q2 .LT. 1.) P13Q2=1.0 

2 W1=141RX(PT2/14.696)XQRTH2 
CXX*¥HFAN STALL MARGIN CALCULATIONS 

14STL=FUN2C8,5,5,XBETA2,YPCNLR,Z14STL, BETA F, PCNLR) 
PRSTL=FUN2(9,5,5,XBETA2,YPCNLR,ZPRSTL,BETAF, PCNLR) 
SM=100. X((PRSTLX141R/P13Q2/(14STLX. 0 929X17. 18/. 4536 ) )-l. ) 

PT13=P13Q2XPT2 

TT13=T13Q2XTT2 

RTT13=SQRT(TT13) 

CXXXXXFAN DUCT CALCULATIONS 
PT18= . 99XPT13 
14A18-I41 
PT25=PT2 


,T11 
, PT25 
»WG4 
, PT8 
, V18 


# TT2 
, TT25 
,WFM 
, PT42 
, FGROSS 


,UFACC, FNZ18 . 

, HPX ,BETAFD,- 
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TT25=TT2 

QRTH25=1 . 0/SQRKTT25K1 . 928E-3) 

CKXXXXCOMPRESSOR CALCULATIONS 

PCNHR=NH/154.5/CSQRT(TT25/605.4)) 

W25R=1.013*FUN1C33,12,XPCNH1,ZW25R,PCNHR) 

ETAC=FUN1C34,16,XPCNH2,ZETAC,PCNHR) 

W25=W25RX.068046*PT25*QRTH25 

34 W3= . 9036XW25 
WFM=WFH/3600. 

WG4=WFM+W3+.025*W25 

FREQ2=0.0104XW3 

TAU1=0.843/FREQ2 

FREQ4=0.0240*WG4 

TAU3=0.868/FREQ4 

CXXXXXITERATE FOR COMPRESSOR DISCHARGE STATIC PRESSURE 
CxxxxxxxHEAT SOAK DYNAMICS 
KPS3=0 

35 ITR3S=ITR3S+1 
PS3=PS3GS 

IFCPS3.LT. 1.0) PS3=1.0 
PT3 = 1 . 057XPS3 
P3QP25=PT3/PT25 

T3Q25=1 .+CP3QP25XX. 2857-1. J/ETAC 

TT3P=T3Q25*TT25 

T3LDT=(TT3P-T3L)XFREQ2 

TT3=TAU1*T3LDT+T3L 

TT4P=53233.X(WFM/W3)+159.+TT3 

T4LDT=(TT4P-T4L)XFREQ4 

TT4=TAU3XT4LDT+T4L 

PS3=.0812XWG4XSQRT(TT4) 

ERR=PS3-PS3GS 

CALL CNVRG1C8, ERR, PS3,PS3GS, .0002, ITR3S,KPS3) 
IFCKPS3 . EQ . 1 )GO TO 40 
GO TO 35 
40 CONTINUE 
PT4=PS3 

CxxxxxHIGH AND LOW TURBINE EXIT CONDITIONS 
TT42 = . 7897XTT 4 
WG8=W25+WFM 
PT42= . 30 OXPT 4 
PT8 = . 01533XWG8+P0 

TT8Q42=1.-(1.-CPT8/PT42)XX.248)X.866 
TT8=TT42*TT8Q42 
CXXXXXFAN ROTOR DYNAMICS 

HP1T=.3817XWG8X(TT42-TT8) 

HPF=.3623HWA18X(TT13-TT2) 

CXXXXXCORE ROTOR DYNAMICS 

HP2T=.417*WG4X(TT4-TT42) 

HP2C=.369*W3*CTT3-TT25) 

NHDT =37710. /NHx(HP2T-HP2C) 

CxxxxxTHRUST CALCULATIONS 
P0Q18=P0/PT18 

FNZ18=FUN1(35,17,XF18,ZF18,P0Q18) 
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CXXXKXOR SUBSTITUTE FNZ18=S«RT(l.-P0Q18HH(2./7.)) 
IFCFNZ18 .LT. 0.) FNZ18=0. 

V18=109.61*RTT13XFNZ18MCV18 

FGROSS=.0311*WA18*V18 
54 FRAM=W1*V0*0. 031085 
FN=FGROSS-FRAM 

CxxxxxCORE INLET AIR TEMPERATURE SENSOR 

TAUT25=FUN1( 36 , 10 ,XF81 ,ZF81 , 1 . 38«W25 ) 
TT25DT=(TT25-T25SN)/TAUT25 
TH25SN=T25SN/605,4 
RTH25S=SQRT(TH25SN) 

CXKKXXFUEL CONTROL 

NHD=154.5*PCNHRDXRTH25S 
NHD=AMIN1(NHD, 17000.) 

NHD=AMAX1(NHD,8795.) 

NHERR=NH-NHD 

CONST=SQRTCWFH)MNH*1.494E-6 
WFPR0P=-2.29*NHERR*C0NST 
C0N=-3 21 

IF(WFC2.NE.WFGV) CON=0. 

WFIDT=NHERRXCON*CONST 

MFGOV=WFPROP+UFI 

WFSCH=FUN2< 37 , 24 , 5 ,XNH , YTT25 , ZUFSCH , Nil , Y iei/SH ) 
WFACMX=-.0134XNH+234.1 
IF(UFSCH.GT.WFACMX) WFSCH=WFACMX 
IF(PS3. LE.200 . ) PC=PS3 

IF(PS3. GT. 200 . ) PC=200.+(PS3-200. )«(~2. ) 
WFACC=UFSCHXPC 
WFDEC= . 385XWFACC 
WFC1=AMIN1<WFACC* WFGOV) 

WFC1=AMAX1(WFC1,UFDEC) 

WFC1=AMIN1(WFC1,4438.) 

UFC1 =AMAX1 ( WFC1 ,295.858) 

WFC2=WFC1 

WFGV=WFGOV 

WFHDT=(WFCl-WFH)/.05 
CXXXXXFAN PITCH ACTUATOR 

BETADT=10 .X(BETAFD-BETAF) 

IFCABS(BETADT) . GT . 100 . ) BETADT=SIGN<100 . , BETADT) 

RETURN 

END 
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Subroutine CNVRG1 (N, ERR , XOLD , XNEW , TOL, ITR, IND) 


A convergence subroutine called from the Iteration loops of the ENGNyD 

subroutines to find compressor-exit static pressure (PS3). The 

agruments are as follows: 

N ■ parameter used to keep track of what variable is being iterated 
on when CNVRG1 is used for more than one iteration loop. 

ERR ■ difference between latest calculated value and latest guessed 
value. 

XOLD - latest calculated value (in ENGNyD) 

XNEW ■ latest guessed value (from CNVRG1) 

TOL ™ sets tolerance on largest difference between latest guessed 
value and latest calculated value within which convergence is 
assumed to occur. 

IND » Index set in CNVRG1 indicating whether or not convergence is 
complete. 

ITR ■ number of iteration loops (maximum of 50,000 per pass through 
ENGNyD) . 
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XHCNVRG1, 09/19/80 12:38:17 


100 


SUBROUTINE CNVRG1 CN, ERR,XOLD,XNEW, TOL , ITR, IND) 

200 


DIMENSION XV1C10),XV2C10),XR1(10),XR2C10> 

300 


IFdTR.GT. 50000) CALL EXIT 

400 


RELTOL=TOL*ABS(XOLD) 

500 


IF(ABSCERR).LT.RELTOL) GO TO 30 

600 


XV2CN)=XV1CN> 

700 


XR2(N)=XR1(N) 

800 


XV1(N)=X0LD 

900 


XR1 (N)=ERR 

1000 


IFCITR.GE.2) GO TO 20 

1100 


XNEU=XNEW*1.04 

1200 


RETURN 

1300 

20 

IF(XRKN). NE.XR2(N) ) GO TO 25 

1400 


XNEW=XNEW*1.01 

1500 


RETURN 

1600 

25 

XNEW=XV1CN)-XR1(N)*(XV2(N)-XV1(N))/(XR2(N>-XR1(N>> 

1700 


RETURN 

1800 

30 

IND=1 

1900 


XNEU=XOLD 

2000 


RETURN 

2100 


END 
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Block Data 


BLOCK DATA contains some miscellaneous constants and coefficients as 
well as all of the tables of data that are functions of one or two 
variables. These data Include items such as: ambient temperature and 
pressure as functions of altitude; fan airflow and temperature and 
pressure ratios as functions of fan speed and fan-blade pitch angle; 
and core-compressor airflow and efficiency as functions of core speed. 
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BLOCK DATA 
REAL MNO 

COMMON /XXI/ ALT,MN0,AE8,CV8,CV18,PS3GS>PT8GS,QX»ICF 
COMMON /XX2/ XALT(8) »ZPFA(8 ) ,ZTFA(8 ) , XPCNH2( 16 ) , ZETACC 16 ) , - 
AXPCNH1 ( 12 ) , ZW25RC 12) >XF18 (17)»ZF18(17) , XF81C10 ) , ZF81 (10),YY(50), DYDT(50 ) 
COMMON /FAN/ XBETAL7 ) , YPCNLRC5 ) ,ZFFLOWC 7 , 5 ) ,ZFTPR(7,5> ,ZFTTR<7 5 ) 

COMMON /HUB/ XBETA1L8) ,ZFHPR(8 , 5) ,ZFHTRC8 , 5 > ' 

COMMON /CORE/ XNH<24> ,YTT25<5) ,ZWFSCHC24,5) 

COMMON/ FA NS TL/XBETA2(5) >ZWSTL(5>5) *ZPRSTLC5»5) 

DATA ALT,MN0,AE8,CV8,CV18/0.0,0.0,371.7, .997, .985/ 

DATA PS3GS, PT8GS/196 .8, 16 .31/ 

DATA XF18/.5, .55, .6 , .65, . 7, . 75, .8, .85, . 9, . 92, . 94, . 95, . 96, . 97, . 98, . 99, 1 ./ 
DATA ZF18/.4239, .3963, .3685, .3403,. 3113,. 2809, .2485, .2130, .1722,- 
A . 1534, . 1324, .1206, .1077,. 0 931,. 0759, .0536,0.0/ 

DATA XF81/ 9. 5, 11. 4, 15. 2, 19. 0,22. 8, 28. 5, 38. 0,47. 5, 57. 0,76.0/ 

DATA ZF81/3. 06, 2. 72, 2. 22, 1.89, 1.64, 1.38, 1.13, 0.97, 0.85, 0.72/ 

DATA XPCNH1/65. ,70. ,75. ,80. ,82. 5, 85. ,90. ,92. 5, 96. ,100. ,105. ,107. 5/ 

DATA ZW25R/12. 2, 14. 0,16. 2, 18. 9, 20. 8, 24. 2, 29. 7, 32. 9, 37. 2, 41.2,- 
A 44.1,44.9/ 

DATA XPCNH2/65. ,70. , 78 . 5 ,82 . 1 ,85 . ,86 . 69,87 . 28, 90 . 07, 91 . 13, 92 . 78, - 
A 94.39,95.28,97.65,98.35,102.23,106.39/ ' 

, DA I^ ZETAC/ * 78 ' * 79 ' * 8, •«!* -82, .83, .832, .839, .842, .844, .845, .846,- 
A .844, .843, .833, .804/ 

DATA XALT/O . 0 , 5 . E3> 1 . E4 > 1 . 5E4 , 2 . E4 , 2 . 5E4 > 3 . E4 > 3 . 5E4/ 

DATA ZPFA/14. 696, 12. 228, 10. 107, 8. 294, 6. 754, 5. 453, 4. 365, 3. 458/ 

DAT A 2TFA/518. 67, 500. 85, 483. 03, 465. 21, 447. 39, 429. 57, 411. 75, 393. 93/ 

DATA XBETA/-20. ,—15. , — 10. ,— 5. ,0. ,5. ,7.3/ 

DATA YPCNLR/70. ,80. ,90 .,100. ,110./ 

DATA ZFFLOW/- 

A 295. ,342. ,389. ,435. ,476 . ,504 . . 513 . ,- 
B 336. ,389. ,442. ,495., 543. ,577. ,589. ,- 
C 378. ,437. ,4 97., 556. ,607. ,650. ,666.,- 
D 421. ,485. ,549. ,612. ,669. ,719. ,739. ,- 
E 463. ,530. ,597., 664., 716. ,750. ,761./ 

DATA ZFTPR/- 


A 

1.039, 

1.049, 

1.063, 

1.080, 

1.097, 

1.110, 

1.114, 

B 

1.049, 

1.064, 

1.083, 

1.105, 

1.127, 

1.145, 

1.152, 

C 

1.060 , 

1.081, 

1.105, 

1.133, 

1.160, 

1.187, 

1.200, 

D 

1.077 , 

1.101, 

1.130, 

1.164, 

1.199, 

1.232, 

1.246, 

E 

1.092, 

1.122, 

1.156, 

1.196, 

1.231, 

1.255, 

1.262/ 

DATA ZFTTR/- 







A 1.0135, 1.0172, 1.0211, 1.0253, 1.0300, 1.0352, 1.0381, 
B 1.0175, 1.0223, 1.0273, 1.0328, 1.0390, 1.0463, 1.0503, 


1.0215, 1.0281, 1.0347, 1.0416, 1.0490, 1.0581, 1.0633, 
1.0266, 1.0345, 1.0427, 1.0513, 1.0605, 1.0714, 1.0775, 


1.077, - 


DATA XBETA1/-20 
DATA ZFHPR/- 

. ,-13. 

, -9 . , -7 

. ,-4.,l 

. ,4. ,7. 

3/ 

1.059, 

1.065, 

1.070, 

1.074, 

1.078, 

1.081, 

1.081 

1.086 , 

1.089, 

1.092, 

1.094, 

1.100, 

1.105, 

1.105 

1.109, 

1.115, 

1.119, 

1.122, 

1.128, 

1.136, 

1.141 

1.140, 

1.141, 

1.144, 

1.151, 

1.162, 

1.175, 

1.187 

1.167, 

1.175, 

1.177, 

1.180, 

1.200 , 

1.213, 

1.224 
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DATA ,09/19/80 12:39:12 


5300 

5400 

5500 

5600 

5700 

5800 

5900 

6000 

6100 

6200 

6300 

6400 

6500 

6600 

6700 

6800 

6900 

7000 

7100 

7200 

7220 

7225 

7230 

7235 

7240 

7245 

7250 

7255 

7260 

7265 

7270 

7275 

7280 

7300 


DATA ZFHTR/- 
A 1.0214, 1.0241 
B 1.0294, 1.0325 
C 1.0374, 1.0409 
D 1.0469, 1.0504 
E 1.0570, 1.0610 
DATA XNH/9000 . , 
A 12500. ,13000., 
B 16200. ,16250., 
DATA YTT25/394 . 
DATA ZWFSCH/18. 
A 26.7,27.5,27.9 
B 13.15,6. ,18.5, 
C 25.4,25.7,25.9 
D 13.15,6. ,18.2, 
E 23.3,24.1,24.7 
F 13. 15, 6. ,17. 8, 
G 23.25,23.6,23. 
H 13. ,6. ,17.4,16 
I 22.4,22.9,23., 
DATA XBETA2/-17 
DATA ZWSTL/ - 


L, 1.0264, 1.0282, 1.0309, 1.0342, 1.0360, 1.0376, 

>, 1.0352, 1.0369, 1.0395, 1.0445, 1.0479, 1.0516, 

>, 1.0443, 1.0463, 1.0498, 1.0561, 1.0605, 1.0653, 

►, 1.0536, 1.0563, 1.0621, 1.0702, 1.0761, 1.0811, 

I, 1.0642, 1.0680, 1.0755, 1.0844, 1.0906, 1.0994/ 

9500.. 10000..10500..10750..11000..11500..12000.,“ 

13 500. . 13750. .14000. .14200. >15800. ,16000. ,16100.,- 

16300.. 16400.. 16500.. 17000.. 17500./ 
67,459.67,518.67,584.67,624.67/ 
8,18.2,17.1,16.4,16.25,16.35,17.7,19.5,21.7,24.5,- 
{ ?8 • » 28 . ,27 .6,26. ,24.4,23.6,22.94,21.62,20.3,- 

17.7.16.7.16.. 15.9.16..17..18.7.20.8.22.7.24.6,- 
',25.9,25.8,25.6,24.275,23.61,22.95,21.62,20.3,- 
17.2,16.25,15.65,15.6,15.65,16.5,17.9,19.9,21.8,- 

5.25. . 25. .25. .24.8.24.6.23.88.23.17.21.73.20.3,- 

7,23.7,23.7,23.7,23.55,23.475,23.4,21.85,20.3,- 
.4, 15. 6, 15. 1,15. ,15. 1,15. 75, 16. 9,18. 4, 20., 21. 7»- 

23. . 23. .23. .23. .22. 9. 22. 8. 22.. 20. 3. 13. 15. 6./ 

.2, -4. 1,0. ,4. ,7.3/ 


A 72.3,96.8,101.6,116. ,132.8, - 
B 83.7,111.2,117.2,137.2,154.4, - 
C 96.7,127.6,138. ,156. ,173.2, - 
D 110. ,145. 6, 156. 8, 172. 4, 192. 4, - 
E 123.3,165.6,174.8,190.4,211.5/ 
DATA ZPRSTl/ - 


A 1.0883,1.1120, 1.1124,1.1128,1.1128, - 
B 1.1153,1.1448,1.1432,1.1464,1.1456, - 
C 1.1483,1.1888,1.1852,1.1828,1.1880, - 
D 1.1840,1.2240,1.2408,1.2352,1.2536, - 
E 1.2233,1.3040,1.3024,1.3072,1.3200/ 

END 
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FUNCTION FUN1 (N,NXP,XX.ZZ,XIN) 

Function routine used in DVTOL and ENGNyD routines to linearly 
interpolate tables of data having one dependent and one independent 
variable. Results are obtained by extrapolation if the range of the 
table is exceeded. . The arguments are as follows: 

N - a parameter used to keep track of which data table is being 
interpolated. 

NXP ■ number of XX, ZZ pairs in the table. 

XX ■ the table of independent variables. 

ZZ ■ the table of dependent variables. 

XIN - the value of the independent variable for which the dependent 
variable is desired. 
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100 

200 

300 

400 

500 

600 

700 

800 

900 

1000 

1100 

1200 

1300 

1400 

1500 

1600 

1700 

1800 

1900 

2000 

2100 

2200 

2250 

2300 

2400 

2500 

2600 

2700 

2800 

2900 


FUNCTION FUN1(N,NXP,XX,ZZ,XIN) 

COMMON /FMEMR/IX( 60 ) , JY(60 ) , IERRC60 ) 

DIMENSION XX(NXP) ,ZZ(NXP) 

CXXXXXTEST FOR X IN PREVIOUS INTERVAL 
I=IX(N> 

IF(XIN-XX( I) > 120 » 200,110 
110 IF(XIN-XX(I+1 1)200,140,140 
CxxxxxcoUNT DOWNxxxxx 

120 IF(XIN~XX(1 ) >160, 160 » 130 
130 1=1-1 

IF(XIN-XXd) >130,200, 200 
CxxxxxcoUNT UPxxxxx 

140 IFCXIN-XXCNXP) >150 ,170,170 
150 1=1+1 

IFCXIN-XXCI+1 >>200,200,150 
160 1=1 

GO TO 180 
170 I=NXP-1 

180 IFC IERR(N> >200 , 190,190 
190 WRlTE(6,400)N,XIN 
IERR(N>=-1 

CXXXXXINTERPOLATE FOR ANSWERXXXXX 
100 CONTINUE 

200 XFRAC=(XIN-XX( I> )/(XX( 1+1 )-XX( I ) ) 

FUN1=ZZ( I >+XFRACX(ZZ( 1 + 1 )-ZZ( I > > 

IX(N> =1 
RETURN 

400 FORMAT (1H0,12H FUNCTION NO.,I3,20H INPUT OUT OF RANGE, 
12X, 6HXIN = ,G12.4> 

END 
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FUNCTION FUN2 (N,NXP.NYC t XX,YY. ZZ.XIN.YIN) 


A function routine used in DVTOL and ENGNyD subroutines to linearly 

Interpolate tables of data having one dependent and two independent 

variables. Results are obtained by extrapolation if the range of the 
table is exceeded. The arguments are as follows: 

N " a parameter used to keep track of which data table is being 

interpolated. 

NXP = number of XX,ZZ pairs in the table of a constant value of YY. 
NYC = number of YY's. 

XX «* the table containing one set of independent variable. 

YY “ the table containing the other set of independent variable. 

ZZ " the table of dependent variables. 

XIN,YIN = values of the two independent variables for which the 
dependent variable is desired 
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100 

200 

300 

*00 

500 

600 

700 

800 

900 

1000 

1100 

1200 

1300 

1*00 

1500 

1600 

1700 

1800 

1900 

2000 

2100 

2200 

2300 

2*00 

2500 

2600 

2700 

2800 

2900 

3000 

3100 

3200 

3300 

3*00 

3500 

3600 

3700 

3800 

3900 

*000 

*100 

*200 

*300 

**00 

*500 

*600 

*700 

*800 

*900 

5000 

5100 


FUNCTION FUN2(N»NXP,NYC,XX, YY,ZZ,XIN»YIN) 

COMMON /FMEMR/IX(60 ) , JYC60 ) , I ERR C 6 0 > 

DIMENSION XX(NXP),YY(NYC), ZZ(NXP, NYC) 

I = IXCN) 

J = JYCN) 

CXXXXXTEST FOR X IN PREVIOUS INTERVALXXXXX 
IF(XIN-XXCI)) 120,200,110 
110 IF(XIN-XX( 1+1 ) ) 200,1*0,1*0 
CXXXXKCOUNT DOWNXXXXX 

120 IF(XIN-XXCll) 160,160,130 
130 I = 1-1 

IF(XIN-XX(I)) 130,200,200 
CXXXXKCOUNT UPXXXXX 

1*0 IFCXIN-XX(NXP) ) 150,170,170 
150 1=1+1 

IF(XIN-XX(I+1)) 200,200,150 
160 I = 1 

GO TO 180 
170 I = NXP-1 

180 IF( IERR(N) ) 200,190,190 
190 WRITE(6 , *00 ) N,XIN,YIN 
IERR(N) = -1 

CxxxxxTEST FOR Y IN PREVIOUS INTERVALXXXXX 
200 IFCYIN-YY(J)) 220,300,210 
210 IF(YIN-YY( J+l ) ) 300,2*0,2*0 
CXXXXKCOUNT DOWNXXXXX 

220 IFCYIN-YYC1)) 260,260,230 
230 J = J-l 

IF( YIN-YYC J ) ) 230,300,300 
CXXXXKCOUNT UPXXXXX 

2*0 IFCYIN-YY(NYC) ) 250,270,270 
250 J = J+l 

IF( YIN-YY( J+l ) ) 300,300,250 
260 J = 1 

GO TO 280 
270 J = NYC-1 

280 IFCIERR(N)) 300,290,290 
290 WRITEC6 , *00 ) N,XIN,YIN 
IERR(N) = -1 

CXXXXXINTERPOLATE FOR ANSWERXXXXX 

300 XFRAC = (XIN-XXC I ) )/(XX( 1+1 )-XX( I ) ) 

P1ZZ = ZZ(I, J)+XFRACX(ZZCI+1, J)-ZZ(I,J)) 

P2ZZ = ZZ(I,J+1)+XFRACX(ZZ(I+1,J+1) -ZZ( I , J + l ) ) 

YFRAC = ( YIN-YY (J))/(YY(J+1) -YY ( J ) ) 

FUN2 = P1ZZ+YFRACX C P2ZZ-P1ZZ) 

I X ( N ) = I 
JYCN) = J 
RETURN 

*00 FORMAT (1H0,12H FUNCTION NO.,I3,20H INPUTS OUT OF RANGE, 
12X, 6HXIN = ,G12.*,2X,6HYIN = ,G12.*) 

END 



SUBROUTINE EULERM (DYDT , YY ) 


This subroutine is called from DVTOL and Is used to find the integrals 
YY of the derivative terms DYDT. A modified Euler method is used as 
outlined in reference 7 . 
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XXEULRM , 0 9/19/80 12:36:55 
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SUBROUTINE EULERMCDYDT ,YY) 

CXXXXXMODIFIED EULER METHOD - REF. MCCRACKEN & DORN P322 
COMMON /INTVAL/DT1,DT2,IC0UNT 
DIMENSION DYDTC50) , YY( 50 ) , YYSV ( 50 ) 

ICOUNT=ICOUNT+l 
IF(ICOUNT) 5,5,20 
5 DO 10 1=1,50 
YYSV ( I ) =YY Cl) 

YY ( I ) =YY ( I )+DT2*DYDT ( I ) 

10 CONTINUE 
GO TO 40 
20 DO 30 1=1,50 

YY(I)=YYSVCI)*DT1*DYDT(I) 

30 CONTINUE 
40 RETURN 
END 
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vectoring 
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Figure l - Sketch of research and technology aircraft showing installation of 
propulsion system. 
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Lift/cruise 
turbofan 
engine 
(unit no. 1) 


Figure 2. - Details of research and technology aircraft propulsion system. 






Ambient Fan inlet Fan duct 

Nozzle 

exit 

For propulsion unit no. 

Conditions at stations 2 and 25 are equal (turboshaft engine inlet) 
temperature at stations 13 and 18 are equal (remote fan exit) 


Figure 3, - Propulsion station numbers for LeRC model. 









(a) Lift/cruise turbofan (units 1 and 2). 

Figure 4. - Computational flow diagram of propulsion system model. 
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(fig. 4(a)) 


(b) Turbo shaft/ remote fan (unit 3). 
Figure 4. - Concluded. 
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Figure 5. - Fuel control details. 
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